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Abstract. Modulo the moment asymptotic expansion, the Cesaro and parametric be- 
haviours of distributions at infinity are equivalent. On the strength of this result, we 
construct the asymptotic analysis for spectral densities, arising from elliptic pseudodif- 
ferential operators. We show how Cesaro developments lead to efficient calculations of 
the expansion coefficients of counting number functionals and Green functions. The 
r-». | bosonic action functional proposed by Chamseddine and Connes can more generally be 

validated as a Cesaro asymptotic development. 
G ■ 

o . 

Ch , 1. Introduction 

<-i_i Most approaches to spectral geometry rely on the asymptotic expansion of the heat 

kernel and Tauberian theorems. In this work, motivated by a string of recent papers by 
Connes, we develop spectral geometry from a more fundamental object. According 
to a deep statement by Connes [10], there is a one-to-one correspondence between 
Riemannian spin geometries and commutative real X-cycles, the dynamics of the latter 
being governed by the spectral properties of its defining Dirac operator. On ordinary 
manifolds, gravity (of the Einstein and the Weyl variety) is the only interaction naturally 
described by the K-cyde [1, 27, 28]. 

That is to say, in noncommutative geometry, existence of gauge fields requires the 
presence of a noncommutative manifold structure, whose "diffeomorphisms" incorpo- 
rate the gauge transformations. Connes' new gauge principle points thus to an intrinsic 
coupling between gravity and the other fundamental interactions. In a recent formula- 
tion [7] , the Yang-Mills action functional is replaced by a "universal" bosonic functional 
of the form: 

B4D]=Tr<p(D 2 ), 
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2 ESTRADA, GRACIA-BONDIA and VARILLY 

with cj) being an "arbitrary" positive function of the Dirac operator D. 

Chamseddine and Connes' work on the universal bosonic functional has two main 
parts. In the first one, they argue that B^ has the following asymptotic development: 

oc 

B^D/k]~Y.^^- 2n a n {D 2 ) as A^oo, (1.1) 

n=0 

where the a n are the coefficients of the heat kernel expansion [19] for D 2 and /o = 
J x(f)(x)dx, fi = f Q (f)(x)dx, J2 = 0(0), /3 = — 0'(O), and so on. Then they proceed 
to compute the development for the -KT-cycle currently [9, 32] associated to the Standard 
Model, indeed obtaining all terms in the bosonic part of the action for the Standard 
Model, plus gravity, plus some new ones. Their approach gives prima facie relations be- 
tween the parameters of the Standard Model, in terms of the cutoff parameter A, falling 
rather wide of the empirical mark. In the second part of their paper, they enterprise to 
improve the situation by use of the renormalization group flow equations [2] . This need 
not concern us here. 

Formula (1.1) can be given a quick derivation, by assuming that <p is a Laplace 
transform. This condition, however, will almost never met in practice. In order to see 
that the asymptotic development of B^ cannot be taken for granted, let us consider, as 
Kastler and coworkers have done [6, 26] the characteristic functions 4>a := X[o,A]- This 
looks harmless enough, giving nothing but A^2(A 2 ), the counting number of eigenval- 
ues of D 2 below the level A 2 . However, it has been known for a long time — see for 
instance [24]-- that there is no asymptotic development for the counting functional 
beyond the first term. Therefore equation (1.1), as it stands, is not applicable to that 
situation. 

One of our aims in this paper is to decrypt the meaning of "arbitrary functional" ; 
a related one is to put on a firm footing the development (1.1). Our contribution turns 
around the Cesaro behaviour of distributions, and its relation with asymptotic analysis. 
Most results are new, or seem ignored in the literature; the paper is written with a 
pedagogical bent. 

The article is organized as follows. Section 2 is the backbone of the paper; there the 
Cesaro behaviour of distributions and Cesaro summability of evaluations are examined. 
The distributional theory of asymptotic expansions [15] is summarized. The latter is 
brought to bear by finding the essential equivalence between the Cesaro behaviour and 
the parametric behaviour of distributions at infinity. Also we prove that a distribution 
satisfies the moment asymptotic expansion iff it belongs to /C', the dual of the space of 
Grossmann-Loupias-Stein operator symbols [20]. These results are new, having been 
obtained very recently by one of us [RE, 12]. We try to enliven this somewhat technical 
section with pertinent examples. 

Next we consider elliptic, positive pseudodifferential operators; let H be one of 
those; the functional calculus for H can be based on the spectral density, formally 
written as 8{\ — H). This is arguably a more basic object than the heat kernel, and 
its study is very rewarding. In Section 3, we show that 5(X — H) is an operator- 
valued distribution in K! . With that in hand, one can proceed to give a meaning to 
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the universal bosonic action for a very wide class of functionals. Following some old 
ideas by Fulling [17], insufficiently exploited up to now, we emphasize that the Cesaro 
behaviour of the spectral density for differential operators is local, i.e., independent of 
the boundary conditions. This is practical for computational purposes, as it sometimes 
allows to replace an operator in question by a more convenient local model. 

In Section 4, we reach the heart of the matter: let d,H(x,y', A) denote the distri- 
butional kernel of 5(\ — H); a formula for dn is given and immediately applied to 
compute the coefficients of its asymptotic expansion on the diagonal, in terms of the 
noncommutative residues [38] of certain powers of H. We hope to have clarified in the 
paper that the identification of the higher Wodzicki terms is essentially a "finite-part" 
calculation. The spectral density is actually a less singular object for operators with 
continuous spectra than for operators with discrete spectra, and all of the above applies 
to operators associated to noncompact manifolds: for that purpose, taking account of 
locality, we work with densities of noncommutative residues throughout. We go on to 
extend Connes' trace theorem [8] to noncompact -ftT-cycles. The case of generalized 
Laplacians is then treated within our procedure. 

In the light of the preceding, the last two sections of the paper are concerned, 
respectively, with the counting number and the heat kernel expansions. The counting 
functional -/Vff(A) is treated mainly by way of example. Then we reexamine the status of 
arbitrary smoothing asymptotic expansions, in particular the Laplace-type expansions 
like the Chamseddine-Connes Ansatz. We point out conditions for the expansions to be 
valid without qualification, and to be valid only in the Cesaro sense. Also we exemplify 
circumstances under which the formal Laplace-type expansion does not say anything 
about the true asymptotic development. The Chamseddine-Connes expansion is derived 
and reinterpreted. 



2. Cesaro computability of distributions 

Besides the standard spaces of test functions and distributions, the space K first intro- 
duced in [20] and its dual KJ play a central role in our considerations. Familiarity with 
the properties of K and KJ and with some of their elements will be convenient. For all 
general matters in distribution theory, we refer to [18]. 

As our interest is mainly in spectral theory, we consider Grossmann-Loupias-Stein 
symbols in one variable, almost exclusively. A smooth function <p of a real variable 
belongs to /C 7 for a real constant 7 if (fi( k \x) = 0(\x\ 1 ~ k ) as \x\ — ► 00, for each k G N. A 
topology for /C 7 is generated by seminorms ||</>||fc, 7 = sup xeM {max(l, (a;^ -7 ) \(fi( k \x)\ }, 
and so /C 7 c — ► /Cy if 7 < 7'. Notice that <f>( k ) G /C 7 -fc if (j> G /C 7 . The space K is the 
inductive limit of the spaces K 1 as 7 — > 00. 

Since every polynomial is in /C, a distribution / G KJ has moments 

ix n := (f(x),x n ), neN 

of all orders; this is an indication that / decays rapidly at infinity in some sense. 

Denote by T>' (T) the space of periodic distributions with zero mean. They con- 
stitute a first class of examples: if / G V (T), then, for n suitably large, the periodic 
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primitive with zero mean f n of / of order n is continuous and defines the evaluation of 
/ at G K, by a convergent integral: 

Note that in this case all the moments are zero. 

The algebra K, is normal (i.e., S is dense in /C) and is a subalgebra of the multi- 
plier algebras Cm, -M of 5, respectively for the ordinary product and the Moyal star 
product [16]. Other properties of /C and /C' will be invoked opportunely. The usefulness 
of K, in phase-space Quantum Mechanics lies in the similitude of behaviour of the or- 
dinary and the Moyal product, when applied to elements of /C. The link between both 
appearances of K is still mysterious to us. 

The natural method of studying generalized functions at infinity is by considering 
the parametric behaviour. The moment asymptotic expansion of a distribution [15] 
is given by 

f(\x)~^ -^j^i asA^oo. (2.1) 

fc=0 

The interpretation of this formula is in the distributional sense, to wit 

(f(Xx)A(x)) = f2^^ + 0( J ^) as A^oo, 

k=0 

for each <p in an appropriate space of test functions. Such an expansion holds only 
for distributions that decay rapidly at infinity, in a sense soon to be made completely 
precise; it certainly does not hold for all tempered distributions, as their moments do 
not generally exist. Distributions endowed with moment asymptotic expansions are said 
to be "distributionally small at infinity" . We are not happy with this terminology and 
invite suggestions to improve it. 

On the other hand, the classical analysis [23] notion of Cesaro or Riesz means 
of series and integrals admits a generalization to the theory of distributions, that we 
intend to exploit in this paper. It turns out that Cesaro limits and "distributional" ones 
are essentially equivalent; this will enable us to apply the simpler ideas of parametric 
analysis to complicated averaging schemes. 

We begin now in earnest by introducing the basic concept of Cesaro behaviour of the 
distributions; justification will follow shortly. Assume / G V(W), /3gE\ { — 1, —2, . . .}. 

Definition 2.1. We say that / is of order x@ at infinity, in the Cesaro sense, and write 



f{x) = 0(xP) (C) 



as x —>■ oo, 



if there exists N G N, a primitive /jv of / of order N and a polynomial p of degree at 
most N — 1, such that /jv is locally integrable for x large and the relation 

f N (x)=p(x)+0(x N+l3 ) asx^oo (2.2) 
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holds in the ordinary sense. 

The relation f(x) = o(ar) (C) is defined similarly. The notation (C,N) can be 
used if one needs to be more specific; if an order relation holds (C, N) for some N, it 
also holds (C, M) for all M > N . The assumption (3 ^ — 1, —2, ... is provisionally made 
in order to avoid dealing with the primitives of x , x~ 2 and such (see Section 6 for the 
general case). If (3 > — 1, the polynomial p is arbitrary and thus irrelevant. We shall 
suppose when needed that our distributions have bounded support, say, on the left. In 
that case, we denote by /[/] the first order primitive of / with support bounded on the 
left. When / is locally integrable, then, 



J — oo 



The notation 

fix) = o(x _ °°) (C) as x — > oo 

will mean fix) = 0{x@) (C) for every (3. 

For the proof of the following workhorse proposition we refer to [12]. 

Lemma 2.1. (a) Let f £T>' such that 

f(x)=0(x 13 ) (C,N) as x^oo. 
Then for k = 1, 2, 3, ... we have: 

fW(x) = 0(x f3 - k ) iC,N + k) as x^oo. 
(b) Let f E V such that 

f{x) = Oix 13 ) (C) as x — > oo, 
and let aeK. Provided that a + ft is not a negative integer, we have: 

x a fix) = Oix a+(3 ) (C) as x -> oo. D 

Definition 2.2. We write limj;^^ f(x) = L (C) when f(x) = L + o(l) (C) as x — > oo. 
That is, lim^^oo f[x) = L (C, k) when fk(x) kl/x h = L + oil), for fk a primitive of 
order k of /. 

For example, if / is periodic with zero mean value, there exists n G N and a con- 
)us (thu 
then clearly 



tinuous (thus bounded) periodic function f n with zero mean value such that /A = /; 



fix) — oix °°) (C) as a; — > oo, 
a fact that yields, for / periodic with mean value ao'. 

lim fix) = a (C). 
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Let / G V be a distribution with support bounded on the left and let be a 
smooth function. The following is a key concept of the theory. 

Definition 2.3. We say that the (f(x), <f>(x)) has the value L in the Cesaro sense, and 
write 

(f(x),(f>(x)) = L (C) 

if there is a primitive I[g] for the distribution g(x) = f(x)(f)(x), satisfying 

lim I[g](x) = L (C) as x — ► oo. 

x — >oo 

A similar definition applies when / has support bounded on the right. If / is an 
arbitrary distribution, let / = /i + fi be a decomposition of /, where f\ has support 
bounded on the left and fi has support bounded on the right. Then we say that 
(f{x), <j>{x)) = L (C) if both (fi(x), <j>{x)) = L % (C) exist for % = 1, 2 and L = L x + L 2 : 
this definition is seen to be independent of the decomposition. 

For instance, let / be a periodic distribution of zero mean and let fi, f2, ■ ■ ■ , fn+i 
denote the periodic primitives with zero mean of /, up to the order n + 1. Then 

x n h(x) - nx^fcix) + n(n - l)x n ~ 2 f 3 (x) -■■■ + (-l) n n! f n+1 (x) 

is a first order primitive of x n f(x), and since fi(x) = o(x~°°) (C) for % — 1, . . . ,n as 
x — > oo, it follows that 

(f(x),x k ) = (C) 

for all fceN. 

To perceive the point of our hitherto abstract definitions, it is worthwhile to recall 
here briefly the classical theory [23]. Let {a n }^ =1 be a sequence of real or complex 
numbers. Often it has no limit, but the sequence of averages Hn := (oi + • • • + a n )/n 
does. Then people write 

lim a n = L (C, 1). 

n — >oo 

If Hn still does not have a limit, then one may apply the averaging procedure 
again and again, hoping that eventually a limit will be obtained. There are two main 
procedures to perform such higher order averages: the Holder means and the Cesaro 
means. The Holder means are single-mindedly constructed as 

zr(k-l) rr(fc-l) 

H (k) ._ U l ^ <- tin 



II 



and lim^^oo H n = L is written 



lim a n = L (H, k). 
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The properly named Cesaro means are denned as follows: let A n :— a n and define 
recursively A n = A\ + ■ ■ ■ + A n . If lim n ^ 00 k\ A n /n k = L, we write 

lim a n = L (C, k), 

n — >oo 

so that the (C, 1) and the (H, 1) notions are identical. The Cesaro limits have nicer ana- 
lytical properties. The good news, at any rate, is that both procedures are equivalent: 

lim a n = L (C, k) •<=>- lim a n = L (H, k). 

n — >oo ' n — >oo 

One uses the simpler notation lim n ^oo a n = L (C) if lim n ^oo a n = L (C, k) for some 
keN. 

A third averaging procedure is equivalent to Cesaro' s, the so-called Riesz typical 
means. For real n, one writes 

lim a n = L (R, k, n) 

if 

lim - VVl ) a n = L. 

n</Lt 

Riesz originally studied this formula for integral //, but the means have more desirable 
properties with \i real. Now, one may study the summability of a series Y^nLi a n by 
studying the generalized function of a real variable f(x) = J2n°=i an ^( x ~ n )- The 
definition of Cesaro limits of distributions is tailored in such a way that (/, 1) (C) 
and Y^n-i a ™ (^) coincide: a primitive of order k of Y^n-i ttn ^( x ~ n ) * s gi ven by 
fk( x ) — 12 n <x( x ~ n ) k ~ la n/(k — 1)! Note that one could consider distributions of the 
form h(x) = XmLi a ™ ^( x ~ P n )i with p n | oo; this gives rise to the (R, k,p n ) means. 
In summary, we have demonstrated the following equivalence. 

Theorem 2.2. The evaluation 

I oo > 

^2a n 5(x -n),<j)(x) ) = L (C) 

n=l I 

holds iff J2n°=i a ™ ^( n ) = L in the Cesaro sense of the theory of summability of series. 

D 

In the same vein: 

Theorem 2.3. If f is locally integrable and supported in (a, oo), then 

(f{x)^{x)) = L (C) 

if and only if 

/>oo 

f(x)<f)(x) dx = L 
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in the Cesaro sense of the theory of summability of integrals. □ 

As shown below, if / E K,' and (f) G /C, then the evaluation (f(x),<j)(x)) is always 
(C)-summable. We pause an instant to show by the example just how useful is the 
concept of Cesaro computability of evaluations. An interesting periodic distribution is 
the Dirac comb Xl^L-oo ^( x ~~ n )- ^ s Inean value is 1; therefore 



oo 



^ S(x-n) = l + f(x), (2.3) 

n= — oo 

with / G £>o(T). The distributions 



OO oo 

J2S(x-n)-H(x-l), £. 

n=l n=l 



S(x — n) — H(x — 1), > 5(x — n) — H{x), 



oo 


<>oo 

- / 0(x) dx 



where H is the Heaviside function, belong to /C'. In effect, take a function <p\ G /C such 
that 0i (x) = 1 for x > 1/2, 0i(x) = for x < 1/4. Then ^(^(E^-oo 5(x - n) - l) 
only differs from Y^n-i $( x ~ n ) ~ ^{x — 1) or Y^n-i $( x ~ n ) — H(x) by a distribution 
of compact support. 

It follows that the evaluation 

oo 

y S(x — n) — H(x — 1), 4>{x) 

n=l 

is Cesaro summable whenever <f> G /C. Now, x a does not belong to K, unless a G N, 
but the previous argument, using <p a (x) = <fii(x)x a , allows us to conclude that the 
evaluation 

/ oo 

Z(a) := / J2 6 ( x -n)- H(x -l),x c 

is (C)-summable for any a G C. Also, Z(a) is an entire function of a, since <p a is. We 
find a formula for Z(a) by observing that if 3? a < — 1 then the evaluation is given by 
the difference of a series and an integral, so that 

00 POO -1 

Z(a) = V n Q - / x a = ((-a) + , Ka < -1. 

£1 ^ a+1 

We have learned a simple proof that Riemann's zeta function is analytic in C\ {1}, 

with residue at s = 1 equal to 1, and one realizes that the evaluation of the ( function 

can be done by Cesaro means (it is only because the zeta function is the outcome of 

a regularization process that it is useful for renormalization in quantum field theory). 

The evaluation (XmLi ${ x ~ n ) ~ H(x), x a ) is slightly more involved. However, we may 

write [14]: 

00 > „i 

V 8{x -n)- H(x), x a ) := Z(a) - F.p. / x a dx, 

n=l / J ° 
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where F.p. stands for the Hadamard finite part of the integral. Now, 

,1 x 

F.p. / x a dx = , a 7^ —1, 

Jo « + l 



therefore, if a 7^ — 1, 



00 „ c 

C (-a) = J> a - F.p. / 

n=l Jo 



x a dx (C), 



in the sense that 



Lxj 



Vn Q -F.p. / t a dt I (C). 
^1 Jo J 

This formula gives a nice representation for (,{a) when 3? a < 1. For instance, £(0) = 
—1/2 simply because the fractional part {x} = x — [x\ of x is periodic of mean 1/2. 
For cm = — 1: 

( L X J i>x \ 

Vfi- td* I = lim (i|x|(|a;J + 1) - \x 2 ) (C); 
^ io / x ^°° 

we find that 

(a - {s})(s - {x} + 1) x 2 _ {x} 2 - {x} x(l - 2{x}) _ 1 

2 2 2 2 12 l ' y ' 

since (1 — 2{x}) and ({x} 2 — {x} + l/6) are periodic of mean zero; we get £(— 1) = —1/12. 
Also, the logarithm of the "functional determinant" can be obtained by this method: 

( L X J i>x \ 

Vlogn- / Xogtdt) (C), 

on using Lemma 2.1. Stirling's formula gives 

xlogx — x — log( [xj !) = xlogx — x — ([x\ + i) log[xJ + |_ccj — log v27r + 0(x _1 ) 

= -xlog(l-i^i) -{x} + ({x}-i)logLxJ -ilog(27r)+0(x- 1 ) 
= -ilog(27r) + 0(x- 1 ) (C) 

since aTog(l — x _1 {x}) + {x} = 0(x _2 ) and ({x} — |) is periodic of mean zero. From 
this it follows that C'(0) = — \ log(27r). This business of Riemann's zeta function is not 
merely amusing; it will be useful later. 

We make ready for the main equivalence result. 
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Theorem 2.4. Let f e V . If a > -1 then 

f(x)=0(\x\ a ) (C) as x^±oo (2.4) 

if and only if 

f(Xx) =0(X a ) as A^oo (2.5) 

in the topology of V . If — j — 1 > a > —j — 2 for some j G N, then (2.4) holds if and 
only if there are constants /aq, . . . , fij such that 

f(Xx) = ^ ^+i + 0(X ) 

fc=0 

in the topology ofV as X — > oo. 

Proof. We prove the theorem in the case / has support bounded on the left. The general 
case follows by using a decomposition / = /i + /2, where /i has support bounded on 
the left and fi has support bounded on the right. First we have to clarify the meaning 
of (2.5). It is a weak or distributional relation: we write f(x, X) = 0(X a ) as A — > oo 
whenever 

(f(x, A), <f>{x)) = 0(X a ) as A ^ oo, 

for all (j) E T>. Note that this yields 

•df(x,\) 



(x)\ = -(f(x,X),<t>'(x)) = 0(X a ). 



dx 

Now, if (2.5) holds, there exists A^ such that the primitive of order iV of f(Xx), with 
respect to x, exists and is bounded by MA Q , say for |x| < 1 and A > Ao- We have then 
a primitive /at of order iV of f(x), such that 

\f N (Xx)\ < MX a+N , \x\ < 1, A > Ao. 

Taking x = 1 and replacing A by x we obtain 

\f N (x)\<Mx a+N , £>A , 

and thus 

f(x) = 0{x a ) (C, N), asi->oo. 

Reciprocally, assume a > — 1 and f(x) = 0(x a ) (C,N), as x — ^ oo. Then, if f^ is the 
(locally integrable for x large) primitive of order N of f with support bounded on the 
left, an obvious estimate gives /at (Ax) = 0(X a+N ), as A — * oo, and on differentiating 
N times with respect to x one obtains X N f(Xx) = 0(X a+N ), so that (2.5) follows. 

The case when a is nonintegral and less than —1 is more involved, as one has to 
deal with the polynomial p in (2.2). Then one shows that the moments 

(f(x),X k ) = Ll k (C) 

up to a certain order exist, those being essentially the coefficients of p. For the gory 
details, we refer once again to [12]. □ 

A characterization of the distributions that have a moment asymptotic expansion 
follows. 
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Theorem 2.5. Let f e V . Then the following are equivalent: 

(a) /e/C'. 

(b) / satisfies 

f(x)=o(\x\-°°) (C) as x ^±00. 

(c) There exist constants fj,o, /j,i, ^2, . . . such that 

fj, S(x) fJki6'(x) /j 2 S"(x) 
/(Ax) — + -^3 as A -> 00 

in £/ie wea/c sense. 

Proof. It is proven in [15] that the elements of K! satisfy the moment asymptotic ex- 
pansion. For the converse, it is enough, as customary, to consider distributions with 
support bounded on one side. We show that if (b) holds, then / e K' for all 7. From 
the hypothesis it follows that f(x) = 0(x~ 1 ~ 2 ) (C) as x — > 00. Thus, for a certain n, 
the n-th order primitive f n of / with support bounded on one side is locally integrable 
and satisfies f n (x) = p(x) + 0(x~ 1 ~ 2+n ) as x — > 00, where the polynomial p has degree 
at most n — 1. We conjure up a compactly supported continuous function g whose 
moments of order up to n — 1 coincide with those of /. If g n is the primitive of order n 
of g with support bounded on the left, then f n (x) — g n (x) = 0(x~ 1 ~ 2+n ). If G /C 7 _ n , 
the integral j^ifnix) - g n (x))<f>(x) dx converges. Hence /=(/„- # n ) (n) + g e JC^. 
The rest is clear. □ 

We get at once a powerful computational method for duality evaluations. 

Corollary 2.6. If f G KJ and <fi G /C, t/ie evaluation (f(x),<f)(x)) is Cesaro summable. 

Proof. It is enough to check for 0=1. But, according to the previous Theorem, if 
/e/C' then /(x) = o(x _ °°) (C) as a; — >■ 00. By the proof of Theorem 2.4, (f(x), 1) is 
(C)-summable. □ 

Fourier transforms are defined by duality and, in general, if / e S' , we cannot 
make sense of f{u) because the evaluation (e lxu , f(x)) is not defined. However, if cj) e /C 
and «^0, Corollary 2.6 guarantees that the Cesaro-sense evaluation (e lxu , <fi(x)) (C) 
is well defined. Thus 

}(u) = ( e lxu , cj)(x)) (C) when e /C, u # 0. 

It is clear that /C C /C'; this follows also from Proposition 4 of [20]. 

Note as well that the moments of / e /C' are (C)-summable. The converse is true: 

Theorem 2.7. Let f e V . If all the moments (f(x),x n ) = [i n (C) exist for n e N, 
i/ien /e/C'. 

For the easy proof, we refer to [12]. □ 

It is clearly worthwhile to characterize spaces of distributions in terms of their 
Cesaro behaviour. Particularly important is the characterization of tempered distribu- 
tions: 
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Theorem 2.8. Let f G V . Then the following statements are equivalent: 

(a) f is a tempered distribution. 

(b) There exists aeR such that 

f(Xx) = 0(X (X ), as A^oo 

in the weak sense. 

(c) There exists a G IR and k EN such that 

f^\x) = o(\xr- k ) (c) 



as x — > oo. 



Proof. Again, it is enough to consider the case when / has support bounded on one 
side. It is well known that if / G S' then there is a primitive F of some order N of 
slow growth at infinity; it follows that f(x) = 0(|;r| a ) (C). The rest is clear, in view 
of the equivalence theorem 2.4 and the fact that distributional order relations can be 
differentiated at will. □ 

We finish by giving several estimates that we will need later. The first one is just 
a rewording of the properties of the distribution (2.3). 

Lemma 2.9. If g G /C and if f_ g(x)dx is defined, then 



E 

n= — oo 



1 f°° 

g{ne) — - I g(x) dx + o(e°°) as e I 0. □ 

^ J — oo 



By the same token: 
Lemma 2.10. If g E /C(M n ) and if f Rn g(x) dx is defined, then 



^2 g{ke) = e~ n I g(x) dx + o(e°°) as e | 0. 



□ 



Lemma 2.11. If g G /C and if / °° g(x) dx is defined, then 



^^ -i poo -'•- 

Y^ 9 ^ = - / g ^ dx +Yl 

n=l £ J ° n=0 



C(-n)^"'(0) £ „ +o(£OO) (B£la 



77- ! 



Proof. This follows from the zeta function example. □ 

(Results of this type were used to prove some formulas by Ramanujan in [15].) 
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3. Spectral densities 

Let Ti be a concrete Hilbert space, the space of square integrable sections of an Eu- 
clidean vector bundle over a Riemannian manifold M , and let H be an elliptic positive 
selfadjoint pseudo differential operator on Ti, with domain X. We consider the deriva- 
tive, in the distributional sense, of the spectral family of projectors -Eh (A) associated 

to tf: 

dE H {\) 



d H {\) :-- 



dX 



For instance, if H is defined on a compact manifold, and < Ai < A2 < • • • is the 
complete set of its eigenvalues, with orthonormal basis of eigenfunctions Uj, the kernel 
of the spectral family is given by [25]: 

Xj<\ 

and the derivative is 

dH(X) ■■= Yl \ u j)( u j\ 6 ( X - X j)- 

3 

This spectral density is a distribution with values in C(X,Ti). The defining 
properties of E(\): 

/CO poo 

dE(X), H= XdE(X) 

-OO J — OO 

(in the weak sense) become, in the language of the previous section: 

I=(d H (X),l), H=(d H (X),X). 

The spectral density is used to construct the functional calculus for H. Indeed, 
we can define 4>(H) whenever / is a distribution such that the evaluation (^^(A), /(A)) 
makes sense, by 

0(#):=<d H (A),0(A)), 

with domain the subspace of the x G 7i for which the evaluation (^{y \ cIh(X)x),(J)(X)). 
is defined for all y G Ti. 

Especially, one is able to deal with the "zeta operator" : 

H- s :=(d H (X),X- s ), (3.1) 

(for ^ sp H) , the heat operator: 

e~ tH ■= (d H (A),e- tA ), £ >0 (3.2) 

and the unitary group of H, which is just the Fourier transform of the spectral density: 

U H (t):=(d H (X),e- ux ). (3.3) 
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The useful symbolic formula 

d H (X) = 5(X-H) 

recommends itself, and we shall employ it from now on. 

We want to study the asymptotic behaviour of S(X — H). Let X n be the domain of 
H n and let X^ := H^Li Xn- The fact that X^ is dense has, in view of the theory of 
Section 2, momentous consequences. We have 

H n = (5(X-H),X n ) 

in the space C{X 00y r hC). Hence, S(X — H) belongs to the space /C'(R, C(X 00 ,'H))- There- 
fore the moment asymptotic expansion holds: 

5(Xa -H)~^ - [ -^ I as a - oo, 

n=0 

and <5(A — H) vanishes to infinite order at infinity in the Cesaro sense: 

S(\ - H) = o(|A| _0 °) (C) as|A|^oo. 

Of course, the last formula is trivial when H is bounded. 

The space T>(M) of test functions is a subspace of •too. We can then realize the 
spectral density by an associated kernel du(x, y, A), an element of V(W, V'(M x M)). 
Ellipticity actually implies that dn(x, y; A) is smooth in (x, y). The expansion 

d H (x,y;X<7)~J2^^^ , n+1 — as a-^oo 

n=0 

holds in principle in the space D'(R, V(M x M)). We also get 

d w (a;,y;A) = o(|A|- oc ) (C) as |A| -»• oo (3.4) 

in the space T>'{M x M). Equation (3.4) is the mother of all incoherence principles. 
For instance, passing to the primitive with respect to A, for an elliptic operator on a 
compact manifold with eigenfunctions tjj n , n G N, one concludes: 

J2 Mx)ij n (y) = o(|A|-°°) (C) as |A| -. oo, 

A„<A 

for x 7^ y, which is Carleman's incoherence relation [5]. 

It should be clear that the expansions cannot hold pointwise in both variables x 
and y, since we cannot set x = y in the distribution 8 (x — y). In fact, our interest in 
this paper lies in the coincidence limit du (x,x; A), which is not distributionally small. 
However, it is proven in [13] that, away from the diagonal of M x M, the expansions 
are valid in the sense of uniform convergence of all derivatives on compacta. On the 
other hand, if Hi and if 2 are two pseudo differential operators whose difference over an 
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open subset U of M is a smoothing operator, and if di(x, y; A) and dz(x, y; A) are the 
corresponding spectral densities, then [13]: 

di(x, y; aX) — d,2(x, y; aX) + o(a~°°) as a — *■ oo 

in T>'(U x U). Also, it can be shown that 

di(x,y; A) = ^2(2;, y; A) + o(A _ °°) (C) as A — > 00 

uniformly on compacts of U x U, even at the diagonal. 

We exemplify the reported behaviour with the simplest possible examples. Let H 
denote first the Laplacian on the real line. Its spectral density is 

d H (x, y; A) = - — -= cos(v / A(a; - y)) 

27TV A 

and therefore it is clear that dn(x,x; A) is not distributionally small, but rather 

dn(xj x; A) = -j= + o(A _ °°) (C) as A — *■ 00. 

Let H denote now the Laplacian on the circle; the eigenvalues are X n = n 2 , n = 
0, 1,2,.. ., with multiplicity 2 from n = 1 on, with normalized eigenfunctions ip^ (x) = 

(2n)- 1 / 2 e ±inx . Therefore 

d H (x, y; A) = — (d(\) + 2 ^ cos n ( x ~ V) <K A - n2 )) ■ 

~ ' n=l 

Then 

S^\x-y)5^(X) 



-^ (6(Xa) + 2 ]T cosn(x - y) 5(Xa - n 2 )) ~ J] 

n=l 

inX^R,!?'^ 1 x S 1 )), while 



7! (JJ+ 1 

n=l j=0 J 



as cr — > 00 



_l 00 

— f5(A) + 2V cosn(a:-y)5(A-n 2 ) N ) = o(A -00 ) (C) as A -> 00 

27T V ^— ' / 

n=l 



if x and y are fixed, x ^ y. 
On the other hand, 



1 00 

d H (x, x;X) = — (d{\) + 2 J2 6 ( x 



n 2 ' 



n=l 
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does not belong to /C'(R, C°°(M)). For the first time in this paper, but not the last, we 
have to find out what the Cesaro behaviour of a given spectral kernel is. We shall have 
recourse to a variety of tricks. For now, applying Lemma 2.11 to g(x) := </>(x 2 ), for a 
Schwartz function, say, we get: 



oc 



1 



J2 (P(en 2 ) = —= / x~ 1/2 (p(x) dx - §0(0) + o(e°°) as e | 0. 

n=l 2 V £ JO 

It is then clear that 

dn(x, x; A) = -j= + o(A _ °°) (C) as A — > oo. 

and it is also immediately clear that the distributional and Cesaro behaviour of the 
spectral density and its kernel are exactly the same as in the previous example. That 
the manifold be compact or not and the spectrum be discrete or continuous is immaterial 
for that purpose. If we seek a boundary problem for the Laplacian, say on a bounded 
interval of the line, we obtain still the same kind of behaviour (off the boundary, where 
a sharp change takes place). Note also the estimate: 

E l^(*)| 2 ~— (C) as A^oo. 

±;A„<A 

As an aside, we turn before closing this section to the functional calculus formulas 
and compare (3.2) with (3.3). Obviously e~ f ^'' has an extension belonging to /C, so there 
is no difficulty in giving a meaning to the heat operator. Also, as we shall see in Section 6, 
it is comparatively easy to study the asymptotic development of the corresponding Green 
function as t j 0. One of the motivations of the present approach to spectral asymptotics 
is to define a sense for expansions of Schrodinger propagators and the like, that do not 
possess a "true" asymptotic expansion. 

Such an approach can be based in the following idea: Theorem 2.8 points to a 
rough duality between KJ and S' . Let g G <S'(R) and find a so that g(Xx) = 0(A Q ) 
weakly as A — ► oo. For any G «S(R), the function $ defined by 

<f>(x):=(g(\x),<t>(\)) x 

is smooth for x ^ since $(x) = \x\~ 1 (g(X), </>(Ao: _1 ))a, and satisfies 

& n \x) = 0(1x1°"") as |x| ^oo. 

Therefore, if / G K! with ^ supp/, we can define (f(x),g(\x)) as a tempered 
distribution. 

When G supp /, we need to ascertain independently smoothness of $ at the origin. 
It turns out that, for this purpose, it is enough to demand distributional smoothness 
of g, i.e., the existence of the distributional values g( n \0), in the sense of [31], for 
n = 0,1,2,.... Then g(tH) admits a distributional expansion in C(X 00 ,TC) as t | 0. 
This can eventually lead to a proper treatment of some questions in quantum field 
theory. We say no more here and refer instead to the forthcoming [13]. In Section 6 of 
this paper, results will be stated for g belonging to «S(R); for the rest of the paper we 
will venture outside safe territory only in examples. 
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4. The Cesaro asymptotic development of dn{x,x; A) 

In this section we obtain the asymptotic expansion for the coincidence limits of spectral 
density kernels. We are fortified with the results of the previous section, implying that 
the Cesaro behaviour of the spectral density of pseudodifferential operators is a local 
matter. 

Let A be any pseudodifferential operator of order a positive integer d, with complete 
symbol cr(A), on the Riemannian manifold M. To simplify the discussion, we consider 
only operators acting on scalars; the treatment of matrix-valued symbols presents no 
further difficulty. The noncommutative or Wodzicki residue of A is defined by inte- 
grating (the trace of) the partial symbol cr- n (A)(x,£) of order — n over the cosphere 
bundle {(*,£): |£| = 1 }: 



WresA:= / <j- n (A)(x,£) d£dx. 



Here dx denotes the canonical volume element on M. If M is not compact, Wres A may 
not exist, but there always exists the local density of the residue / s „_i cr- n (A)(x, u) du, 
that we denote by wresA(x). 
We recall that 



a(AB) - a(A)a(B) ~ £ K —^d%a{A)d«a{B) 

|a|>0 



(-1)1 



The kernel kA of A is by definition: 

k A (x,y) := (2ir)- n (e^-^,a(A)(x^))^ 

In particular, on the diagonal: 

k A (x,x) := (2ir)- n (l,a(A)(x,Z)) c (4.1) 

In order to figure out the symbol for a spectral density, we start by considering 
(the selfadjoint extension of) an elliptic operator H with constant coefficients. In this 
case a(H n ) = a(H) n and we assert: 

a(6(X-H)) =S(X-a(H)), 

justified by the identities: 



X n 5(X-a(H))d\ = a(H n ), A = 0,1,2,... 

In the general case of nonconstant coefficients, we make the Ansatz that: 

a(8(X -H))~6(\-a(H ))- qi 6'(X- a(H))+ q 2 5"(X-a(H)) -q 3 5'" (X- a(H )) + ••• 

(4.2) 
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in the Cesaro sense. Computation of f X n a(5(X — H)) dX for A = 0, 1, 2, . . . then gives 

91 = 0; q 2 = \{a{H 2 )-a{H) 2 ) ] q 3 = ±(a(tf 3 ) - 3a(H 2 )a(H) + 2a(tf) 3 ), (4.3) 

and so on. This development, it turns out, gives ever lower powers of A in the asymptotic 
expansion of a(S(X — H)). 

We are interested in explicit formulas for the Cesaro asymptotic development of 
the coincidence limit for the kernel of a positive operator H as A — ► oo. From (4.1) 
and (4.2) with p := cr(H), we get 

d H (x, x- A) ~ (27r)- n (l, S(X - p(x, £)) + ? 2 (x, S"(X - p(x, 0) - • • "> c (<?)• 

In polar coordinates on the cotangent fibres, £ = |£|cj with |u;| = 1, this becomes 

(2n)- n f du (ler- 1 , *(A - P (x, ieiw)) + ©(x, ieiw) <*"(a - P (x, ieiw)) -•••>,*• 

•/ |w|=l 

Hence, if we denote by |£|(x, cu; A) the positive solution of the equation p(x, \£\w) = A, 
we need to compute: 

(2) -n f ^ \^\ n -H^^>) + ^(Q2(x,mx,U-,X)um n - 1 (x,U-,X))---- 

' '"' p'(x,|£|(x,u;;A)u;) l " ) 



n-l 



Write: 

p(a ; ,|£| W )~p d (x, W )|£| d +p d _ 1 (x, W )|£| d - 1 +p d _ 2 (x, W )|£| d - 2 ---. 

To solve p(x, \£\u>) = X amounts to a series reversion. 

In order to see how that is done, let us assume for a short while that H is a first- 
order operator with constant coefficients — for instance, the absolute value of the Dirac 
operator on IR n . We then expect 

\t\( w l \ Po ^ ( w-ij 

\i {x,u;X )~— — -A — -p_i^ A H . 

pi(u/) pi(u>) 

Integration over |o>| = 1 gives 

M^,x; A) ~ (27r)- n (a A"" 1 + ai X n ~ 2 + a 2 A n " 3 + • • •) (C), 

where, clearly, ao = wresH~ n . 

To compute ai, a 2 , . . . we can as well assume that the development of p is analytic 

as |£| — * oo. Let ip(z) := z n ~ x /p'(z), so that 



aoX 71 ' 1 + ai (x)X n - 2 + a 2 (x)X n - s + / if>(\£\(x, u; A)) dw. 
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If T is a circle containing |£|(x,u;; A), wound once around oo, we have the Cauchy 
integral: 

i(\c\t \\\ it -i{\\\ 1 I^(z)p'(z)dz 

m\(^\))=^(p w) = -^-% p{z) _ x 

1 I iKrVCC^K 1 f d( 



2mJ r -i ^(KC-^-A) 27rz7 r - 1 C"+ 1 b(C- 1 )-A)' 

Thus cij(x) = L n -i c j( x , w ) du, where 

2 ™ J\s\=e 



1 I so-l-'ds 1 ^ 



(2ttz)2 Jj s|=£ ' 7r-i C" +1 (P(VC) - V«) 

1 I d( I s n ~ j - 1 ds 



(27TZ)2 J v _ x C n+1 J\s\=e S P (l/0 "I 

i r d( 



2m / r _! c n+ MVC) 



n-j 



which is the coefficient of C n in the expansion of p(l/Q J n . Integrating over \uj\ = 1 
yields thus 



a,j = wresH J n , 



so, finally: 



d H (x,x;X)~ —}—( W resH- n \ n - 1 + wresH- n+1 \ n - 2 + ---) (C), 

(2ii) n 

where the densities of Wodzicki residues are constant for a constant-coefficient operator. 
It is amusing that we have arrived at a version of the classical Lagrange-Biirmann 
expansion [29] , with Wodzicki residues in the place of ordinary residues. 

Notice that a n = 0. This is a very simple "vanishing theorem" (see for instance [3]). 

Returning to the general case, if H is a positive pseudo differential operator of 
order d, then A := H 1 ' d is a positive pseudo differential operator of first order. Setting 
|U = X 1 ^, we have 

S(X-H) = S( fl -A) = -^^-= dX(d _ 1)/d . 
and so 

d H {x, x- X) - — I— (a (x)A (n " d) / d + a 1 (x)X^- d - 1 ^ d + a 2 (x)X^- d - 2 ^ d + • • •) (C). 
a {27t) n 

(4-5) 

Clearly, ao = wres H n ' d . Now, the order of q 2 is at most 2d — 1, therefore its higher 

order contribution to this development is in principle to a\; the order of qs is at most 

3d — 2, so it contributes to a 2 at the earliest, and so on. 

Formula (4.5), obtained through fairly elementary manipulations, is the main result 

of this section. To illustrate its power, we show how to reap from it a rich harvest of 

classical results (with a little extra effort) . 
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Corollary 4.1. (Connes' trace theorem) For positive elliptic pseudodifferential ope- 
rators of order —n on a compact n- dimensional manifold, the Dixmier trace and the 
Wodzicki residue are proportional: 

Dtr H = — — — Wres H. 

n(27i) n 

Proof. Let H be of order d = —n in (4.5). We get 

d H (x,x;X)~ — — wres H(x) A" 2 H (C). 

n (2TT) n 

Assume the manifold is compact. We then know that if is a compact operator. Now, 
heuristically the argument goes as follows: iV^(A) ~ —A -2 , ergo Njj(\) ~ A -1 , ergo 
Xi(H) ~ / _1 . A Tauberian argument can be used at this point [37] to ensure that the 
second asymptotic estimate is valid without the Cesaro condition; and then the result 
follows. But this is by no means necessary. One can steal a look at Section 6 and, by 
approaching step functions by elements of S, prove in an elementary way that for any 
given e > there is 1(e) such that 



1(e) " ' 1(e) ' 

where C = n~ 1 (2n)~ n Wres H. □ 

On a noncompact spin manifold, consider now the Dirac operator on the space 
of spinors L (S). The noncommutative integral of |-D| _n does not exist. However, if 
/ a(x) dx is defined, it is computable by a noncommutative integral: 

Theorem 4.2. Let a be an integrable function with respect to the volume form on M . 
Then 

C n / a(x)dx = — ; — — Wres(a|.D|~ n ), 
n J M K ' n{2ix) n ^ I " '' 

where on the right hand side A is seen as a multiplication operator on L 2 (S). The 
constants are C 2 k = (2n)~ k /k\ and C 2 k+i = n~ k ~ 1 /(2k + 1)!! 

Proof. That follows from Theorem 5.3 of [37] if a is a smooth function with compact 
support. For a positive and integrable, use monotone convergence on both sides; the 
general case follows at once. □ 

The former is a small step in the direction of a theory of X-cycles (or "spectral 
triples" , as they are nowadays called) over noncompact manifolds. 

Corollary 4.3. (Weyl's estimate) Let Nh(X) denote the counting function of H , a 
Laplacian on a compact manifold or bounded region M acting on scalar functions. Then 

N H (X) ~ %^ A -/ 2 , 
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where fi n is the surface area of the unit ball in M, n . 

Proof. The same type of arguments as in Corollary 4.1 work. Indeed, this estimate is a 
corollary of it [37] . □ 

Next consider Schrodinger operators —A + V(x), with symbol p(x,£) = |£| 2 + 
V(x). We can take a slightly different tack and solve the equation p(x,£) = A by 

iei = v / (A-n^)) + - 

Corollary 4.4. (The correspondence principle) For Schrodinger operators: 

See [22], for instance, for the reasons for the terminology. 

A word of caution is in order here. The development (4.5) cannot be integrated 
term by term in general. Consider, for instance, the harmonic oscillator hamiltonian 
H = ^(—d 2 /dx 2 +x 2 ) on R: according to the theory developed here, its spectral density 
behaves as 1/vX If ipm n £ N denote the normalized wavefunctions, then indeed, like 
in Fourier series theory, 

£ *fa)~ — 

n+i<A 

is true and can be independently checked. But wresif -1 ' 2 is not integrable over the 
real line, so one cannot conclude that Njj(X) behaves as vA. Actually, as we saw in 
Section 2, £^° =0 <*(*-("+§)) = H(\)+o(\-°°) (C), so N H (X) = \H{\)+o{\-°°) (C). 
Now, Corollary 4.2 applies, so we have 



2 n 2X 

N H (X) / \/2A - x 2 dx = XH(X) 

27r J-V2X 

precisely as it should. (See the discussion in [30].) 

Consider n-dimensional Schrodinger operators with (continuous) homogeneous po- 
tentials V(x) > 0, V(ax) = t a V(x). The previous formula gives 

N H (X) oc X n/2+n/a f V(x)~ n/a dx. 
is™- 1 

and this means that if the cone {x G M. n : V(x) = } is too big, in the counting number 
estimate we are heading for trouble [36]. But the "nonstandard asymptotics" that might 
then intervene do not detract from the validity of the nonintegrated formula (4.5). 

In the remainder of the section, we focus on the computation of spectral densities 
for Laplacians. Nothing essential is won or lost by considering general vector bundles, 
so we work on scalars. The more general Laplacian operator on a Riemannian manifold 
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is (minus) the Laplace-Beltrami operator A plus potential vector and scalar potential 
terms, with symbol 

P(x,0 = -9 lJ (x){^j + {^%{x)ik + 2Ai(x)^) 

+ {A i {x)A j {x) + i(T^(x)A k (x) - diAj(x)))) + V(x) 
=:-g lJ (x)Uj+B\x)^ + C(x). 

Formula (4.5) would seem to give for this case: 

d H (x, x- A) ~ — i— (ao(x)A(- 2 )/ 2 + ai {x)\^-^ 2 + a 2 (x)A^" 4 )/ 2 + • • •) (C). 

2 {Z7T) n 

In fact, it will be seen in a moment that a\ — 03 = • • • = 0. Also we know already 
that ao(x) = WresA _n / 2 = fi n . Our task is to compute the next coefficients; it is a 
rather exhausting one, whose results can be inferred from the extensive work already 
carried out [19] on heat kernel expansions (see Section 6), so we will limit ourselves to 
the computation of a 2 (a;) to illustrate the relative simplicity of our approach. 

Let n > 3. Write a for g l: '(x)uJiUJj, then b for B l (x)iVi and c for C(x). Our method 
calls for solving for the positive root of a|£| 2 + 6|£| + (c — A) = and substituting this in 
|£| n-1 /(2a|£| +b). In diminishing powers of A, we obtain for the latter the development: 

A (n-2)/2 _ (n-l)6 A fa- 3 )/ 2 + Mn-2)b 2 _ (n-2)c \ A ( n _ 4)/2 + \ {4Q) 



2a n / 2 V 2a 1 / 2 V 8a 

One sees that odd-numbered terms in this expansion contain odd powers of to and 
thus give vanishing contributions, after the integration on the cosphere. Also, the 
contribution of the g 2 term in (4.2) will start at order |n — 2 in A, the contribution of 
03 will start at order |n — 3 and so on: the terms in the asymptotic expansion of the 
density kernels of Laplacian operators differ by powers of A, not of y/\, as one would 
expect on general grounds. 

It is convenient now to use geodesic coordinates at each point; this is justified by 
the nature of the result. In these coordinates Y^Axq) = and we have the Taylor 
expansion 

( rp rp _ 1 t-* 

gij(x) ~ Sij + -Rikij(xo) (x - x ) k (x - x ) 1 + 2J d a g(x ) j as x ->■ x , 

|o|>3 

where Rikij denotes the Riemann curvature tensor. Recall that the Ricci tensor is given 
by Rkj := Xw ^klj an< ^ ^ ne sca l ar curvature by R := Ylkj 9 k ^kj- 
From (4.3) one obtains for g 2 (xo,^) 

\ E ^dU-^M^j + B'ixo^d^l^-g^ixMj + B'ix^ + Cix)). (4.7) 

|a|>0 



ON SUMMABILITY OF DISTRIBUTIONS AND SPECTRAL GEOMETRY 23 

Let us take for a moment A^ = 0. Then in geodesic coordinates B i (xo) = and it is not 
hard to see that the only surviving term in (4.7) is equal to ^Rkj(xo)^ k ^ . Also 6 = 
in (4.6). So, in view of (4.4) we are left with two terms at order A*-™ -4 -*' 2 , to wit: 

du (n-2)C(s ) (w _ 4)/2 
2 

that comes from the third term in (4.6), and the first order contribution of 

d ^ -^{q2(x,mx,u;;X)uj)\^- 1 (x,uj;X)) 
„_i j/(x,\^\(x,u;X)u) 

In effect, q 2 contributes here a factor of order A, so the second derivative in the previous 
formula gives rise to a term of order \( n ~ 4: >/ 2 also. To finish the computation, we use 

duo A ij oo' l uj j = —g^An, 
.,..," i- n 

to get 

a 2 (x ) = {U ~l )Qn (\R{x ) - C(x )). (4.8) 

Notice that for a pure Laplace-Beltrami operator, the contribution to a 2 , when com- 
puted in geodesic coordinates, comes exclusively through the q 2 term. 

It remains to convince ourselves that vector potentials give no contribution at this 
stage. On one hand, the c term in (4.6) would contribute now the extra terms 

- (n ~ 2) "" (AM J+ ^). 



On the other, the term in b 2 in the same formula would contribute a term of the form 
\{n — 2)fi n AMj, and in the computation of q 2 there appears now a term (2i/n) djB-i 
that contributes | (n — 2)Q n djA : > and thereby cancels the rest. Therefore (4.8) stands 
also in that case. 

Actually the coefficients of the Cesaro asymptotic expansion of d(x, x; A) are all 
(local densities of) Wodzicki residues for n odd: a 2 k{x) = wres A~ n ' 2+k (x), for k E N. 
For n even we have a 2 k = wres A _n ' 2+fc only as long as —n/2 + k < (the Wodzicki 
residues of nonnegative powers of a differential operator being of course zero); the 
following coefficients for the parametric expansion are, in our terminology of Section 2 
(further explained in the next two sections), not "residues" but "moments". Note that 
for n — 2, the coefficient a 2 is already a "moment" and cannot be computed by a Cesaro 
development. This strikingly different behaviour of the odd-dimensional and the even- 
dimensional cases is concealed in the uniformity of the usual heat kernel method, but it 
reflects itself in the fact that the corresponding zeta functions have an infinite number 
of poles, corresponding to the residues, in the odd-dimensional case; and a finite number 
in the even-dimensional case. One has [38]: 

Res s = n/2 _ k ( H (s) = ±WvesH k - n / 2 , 
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where 



(h(s) - / (d H (x,x;X),X s ) x dx (Ss>0) 
J M 



is the kernel of the zeta operator (3.1). A direct, "elementary" proof of the essential 
identity between Wodzicki residues and residues of the poles of the zeta functions is 
obviously in the cards, but we will not go further afield here. For a nontrivial use of the 
noncommutative residue in zeta function theory, have a look at [11]. 



5. Cesaro developments of counting functions 

We consider here operators on compact manifolds without boundary and look at the 



behaviour of the counting function 



N(X) := J2 1- 



A,<A 

In order to refresh our intuition, we shall follow a deliberately naive approach and 
temporarily forget some of what we learned at the end of last section. Envisage first 
the scalar Laplacian on T 2 with the fiat metric; then the counting function is given by 
the following table: 



A 





1 


2 


4 


5 


8 


9 


10 


13 


16 


17 


18 


20 


25 


26 




N(X+) 


1 


5 


9 


13 


21 


25 


29 


37 


45 


49 


57 


61 


69 


81 


89 





No doubt, N(X) ~ nX is a reasonable first approximation; but it is also plain that the 
remainder undergoes wild oscillations. The precise determination of this remainder is a 
difficult problem, not unlike the problem of determining the next-to-main term in the 
asymptotic development of prime numbers. 

An even simpler and more telling example is provided by the eigenvalues A; of the 
Laplacian on the n-dimensional sphere. They are given by 



Aj — 1(1 + n — 1) with respective multiplicities mi 



l + n 



n 



l + n-2 



II 



(5.1) 



for / G N. For example, if n = 2, the eigenvalues are 1(1 + 1) and the multiplicities are 
(2/ + 1). The leading term is 



N(X) : A n/2 as A -> oo. 



in 



On the other hand, asymptotically: 



N(X + )-N(X~) 



2/ 



n-l 



(ra-1)! 
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and so 



A (1 - n)/2 (iV(A + )-iV(A-)) 



(ra-1)! 

Plainly, we cannot find an asymptotic formula for N(X) with error term o{X^ n ~ 1 '' 2 ) and 
continuous main term. The example is taken from Hormander's work [24, 25]. 

The foregoing is a "Gibbs phenomenon" related to the lack of smoothness of the 
characteristic function. The problem is "solved" if one is prepared to look at the expan- 
sions in the Cesaro sense. The fact that higher order terms in the asymptotic expansion 
of the eigenvalues of the Laplacian were to be understood in an averaged sense was 
pointed out by Brownell [4] many years ago. 

Going back to tori, consider the distribution of nonvanishing eigenvalues {Xi}f^ =1 
of the scalar Laplacian on an n-dimensional torus T n , with the flat metric. The eigen- 
functions {4>i}fZi can be seen as nonzero smooth functions in M. n that satisfy 

A(/>i + XkPi = 

and the periodicity conditions 

01 (xi + 2/citt, ...,x n + 2k n ii) = 4>i(x u . . . , x n ), 

where the girths of the torus are taken to be 2n in all directions. 

Those eigenvalues are given by A& = k\ + ■ ■ ■ + k^ for k = (k±, . . . , k n ) G Z n , with 
corresponding eigenfunctions (f>k(xi, ■ ■ ■ ,x n ) = e lk ' x . Thus the A; are the nonnegative 
integers qi that can be written as a sum of n squares. The multiplicity of each such 
value is the number of integral solutions of the Diophantine equation q\ = k\ + ■ ■ ■ + k\. 
We wish to compute the terms in the parametric and Cesaro developments of N(X) next 
to leading Weyl term (which in fact for this problem goes back to Gauss): 

jV(A) ~ — X n / 2 as A -. oo. 
n 

To do so, we start with the derivative N'(X); this is nothing but (2n) n d(x,x; A), 
but, as advertised, it is more instructive to forget for a while the discussion in Section 4. 
We have: 



Let 4> G £>(K), let a be a large real parameter and set e = 1/cr, so that s J, 0. Then 

(N'(aX),<p(X)) x =e(N'(x),<l ) (eX)) x = eJ2 ^\ k \ 2 ) 

kei n 

= e 1 - n ' 2 [ <f)(\x\ 2 ) dx + o(e°°) 

/•OO 

= ifi n e 1 -^ 2 / r (n " 2)/2 (/)(r) dr + o(e°°). 



o 
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The third equality is just Lemma 2.10. 
Hence, weakly: 

N'(aX) = iO n a" 1+n / 2 A; 1+n/2 + o^" 00 ) as a -> oo, 
and upon integration 



n 



iV( (T A) = ^A" /2 (7 n / 2 +o((7- 00 ) 
n 



as a — > oo. 



Observe that the constant of integration jjlq vanishes, as do all the other moments. 
Then Theorem 2.4 yields: 

jV(A) = ^A n/2 + o(A-°°) (C) asA^oo. 
n 

Hence the error term, although definitely not small in the ordinary sense, is of rapid 
decay in the (C) sense. 

We turn to examine some cases of spheres. The derivative of the counting function 
for S 2 is N'(X) — ^z^ (2/ + 1) S(X — 1(1 + 1)). To deal with this case, we need a heavier 
gun than Lemmata 2.9-2.11. This is provided by: 

Lemma 5.1. Let f G /C'(IR n ), so that it satisfies the moment asymptotic expansion. If 
p is an elliptic polynomial and <ft E S, then 

A — ' to! 

m=0 

Proof. The proof consists in showing that the Taylor expansion 

A </>( m ) (0)p(x 

m=0 



^— ' rn.i 



to! 



holds not only pointwise, but also in the topology of /C(IR n ). D 

Consider now the distribution 

/(A) :^2A + 1) (f>(A - /) - if (A) Y 

that lies in /C'. Notice that 

</(A), ^(t(A 2 + A))> = V(2J + 1) ^(t(/ 2 + /))-/ (2A + l)0(t(A 2 + A)) rfA 

00 poo 

= V(2Z + 1) <p(t(l 2 + 1))- / 0(t/i) d//. 
1=1 Jo 
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From Lemma 5.1 we conclude that, for e S, 

CO 

(N'(\),<Kt\)) = J2( 21 + !) ^W 2 + 0) 

^0 „■_ n ■?' 



r>oo 



(/(A),(A 2 + A)^)(0)^ 
The parametric expansion of N'(X) is thus 



iV'(A/£) ~ #(A) + «J(A)t + f] ( 1) J M 0) (A) tJ+ i as n q 5 

3=0 J ' 

where the "generalized moments" \ij are given by 

00 /-00 

H = (/(A), (A 2 + A)'") = £(2J + 1)(/ 2 + J ' - / (2A + 1)(A 2 + A)' dX (C). 

It follows that AT' (A) ~ if (A) + o(A-°°) (C) as A -► oo. 

In view of our gymnastics with Riemann's zeta function in Section 2, the compu- 
tation of the fjij presents no difficulties. We obtain 

^o = 2C(-1) + C(0) = -% ^2 = 2C(-5) + 4C(-3) 



3 r ,x ' ,v ' 315 

Mi = 2C(-3) + C(-l) = ~, ^3 = 2C(-7) + 9C(-5) + C(-3) 



15 105 

and so on. On integrating, we get 

N(X/t) ~ ^#(A) + ±ff(A) + -U(A) t + ^'(A) t 2 + • • • as t | 0, (5.2) 

t 6 15 315 

and N(X) ~ A#(A) + § H (A) + o(A _0 °) (C). Note that the A°th order term in the 
Cesaro development for N(X) comes from the first moment. The curvature of a sphere 
§ n is given by R = n(n — 1), so the second term in the development is precisely what 
we had expected. 

We look now at the derivative of the counting function for the Laplace-Beltrami 
operator on S 3 . It is slightly simpler to consider the operator 1 — A, for which we have, 

according to (5.1): JV'(A) = E^o^ + l f 6 ( X ~ ( l + !) 2 )- 
Consider the distribution 



/(A) := (A + l) 2 (f>(A - Z) - If (A + 1)\ 
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lying in KJ . We have: 

</(A), <j>{t{\ + l) 2 )) = YJl + V 2 m + l) 2 ) - / (A + 1) 



l (j)(t(\ + l) 2 )d\. 



One sees that the moments all cancel: (/(A), (A + l) 2j ) = ((-2j - 2) = 0, for j e N. 
Therefore we get simply 

1 Z" 00 
(JV'(A),0(tA))~2p72 / <f>(u)yfiidu as tlO, 

and thus in this case we collect just the Weyl term 

\3/2u/n 

jV(A) ~ - — ^Z (C) as A -► 00. (5.3) 

We may reflect now that the counting number for these Laplacians on § 2 , S 3 behave 
in the expected way for even and odd dimensional cases, respectively. For a generalized 
Laplacian which is the square of a Dirac operator the qualitative picture is the same. 
In particular, the Chamseddine-Connes expansion corresponds to n — 4, whereupon 
the counting functional behaves in much the same way as the one for S 2 . Therefore, 
formal application of the Chamseddine-Connes Ansatz to the characteristic function of 
the spectrum, as done in [6, 26] misses the terms involving 5 and its derivatives — whose 
physical meaning, if any, is unclear to us. 



6. Spectral density and the heat kernel 

Now we tackle the issue of the small-t behaviour of the Green functions associated to an 
elliptic pseudo differential operator H. These are the integral kernels of operator- valued 
functions of H, of the form 

G(t, x, y) = (d H (x, y; A), g(t\)) x 

where g, as already advertised, will in this section belong (or can be extended) to the 
Schwartz space S (i.e., we deal with the standard theory as opposed to the framework 
sketched at the end of Section 3). The basic question is whether G(t,x,y) has an 
asymptotic expansion as t J, 0. In effect, we shall see immediately how to obtain from 
the (C) asymptotic expansion for the spectral density an ordinary asymptotic expansion 
for Green functions. 

The emphasis in recent years has been on Abelian type expansions, the so-called 
heat kernel techniques [19]. It is common folklore that Cesaro summability implies 
Abel summability, but not conversely. As we just claimed, one can go from the Cesaro 
expansion to the heat kernel expansion. The reverse implication does not work quite the 
same. If we know the coefficients of the heat kernel expansion and we independently 
know that a Cesaro type expansion for the spectral density exists, we can infer the 
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coefficients of the latter from the former. But it may happen that the formal Abel- 
Laplace type expansion does not say anything about the "true" asymptotic development. 

For instance, if /(A) := sinAe^ A for A > 0, then liniA^oo /(A) (C) does not exist, 
since no primitive of / can have polynomial order in A. Even so, one can show that 
k(t) = (/(A), e~ tx ) has a Laplace expansion k(t) ~ a_it _1 + ao + a±t + • • • as t I 0, 
that is, liniA^oo /(A) = o_i (A). To get an example of a bounded function with this 
behaviour, one uses the fact that / m (A) = sin A 1 /™ obeys lim,\_ ) .oo / m (A) = (C,N) 
only for N > m, together with Baire's theorem, to construct a bounded function /(A) = 
J2 k>1 2~ k /m fc (A) that does not have a Cesaro limit as A — ► oo, but for which /(A) — ► 
in the Abel sense. 

In order to relate our Cesaro asymptotic expansions with heat kernel developments, 
we need to examine expansions of distributions /(A) that may contain nonintegral pow- 
ers of A. Suppose that {ak}k>i is a decreasing sequence of real numbers, not including 
negative integers, and suppose further that / G S' , supported in [0, oo), has the Cesaro 
asymptotic expansion 

/(A) ~ J2 c fc A ° fc + J2 b i X ~ J ( c ) as A ~* °°- 
fe>i j>i 

It follows from Theorem 32 of [15] and from Theorem 2.5 that / has the following 
parametric development: 

f(aX) ~ J>(.A + r + J>Pf((crA)-^(A)) + £ ( " 1) ^C )(A) (6-1) 

fe>l j>1 m>0 

as a — ► oo, where the "generalized moments" \x m are given by 

Li m = (f(x) -J2 c k< k ~J2 b J P?(x~ J H(x)), x m ) (6.2) 

fc>i j>i 

and where Pf denotes a "pseudofunction" [14] obtained by taking the Hadamard finite 
part, that is: (Pi(h(x)),g(x)) := F.p. J °° h(x)g(x) dx if supp/i C [0, oo). In particular, 



r 

(Pf(x-^(x)),^(x)> = F.p. / 

./0 






izi„(fc)rm ,\ ez, 2 y (fc)(o) (6 - 3) 



h X> Jo *> V £"n ^ ! / £"n 



fe!0" — A; — 1) 

fc=0 fc=o u y 

Notice that taking the finite part involves dropping a logarithmic term proportional 
to (7^ _1 ' ) (0). This has the consequence that Pf(x _J i/(x)) fails to be homogeneous of 
degree —j by a logarithmic term; indeed, 

Pf((aA)"^(aA)) = ^ Pf A^tf A + l J °. ^ J ? - ■ 
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Consequently, 

X^g(X)dX 

+ E /W^(0) t „ (64) 



m>0 



The heat kernel development may be recovered by taking g(X) = e x for A > 0. In 
that case, da k is integral, F.p. J °° X ak g(X) d\ = T(a k + 1) and # (j_1) (0) = (-l) (i_1) . 
From this it is clear that the heat kernel of a pseudodifferential operator may generally 
contain logarithmic terms. Indeed, by harking back to (4.5), on using (6.4) we prove: 

Corollary 6.1. The general form of the (coincidence limit of) the heat kernel for an 
elliptic pseudodifferential operator of order d on a compact manifold M of dimension n 
is given by 

oo 

K(t,x,x)~ Yl l3-n{x)t {j - n)/d + Yl Pj-n(x)t (j - n V d logt + J2r m (x)t m 

j-ngdN + j-nGdN + r=l 



qs t I 0, where 



, , T((n — j)/d) , . 



and similarly for the other coefficients. 

(See [21, Cor. 4.2.7].) 

Now suppose we know a priori that /(A) has a Cesaro asymptotic expansion in 
falling powers of A, and that we also know that $(£) := (/(A), e _tA ) A has an asymptotic 
expansion as t I without logt terms. Then it follows that all bj = in (6.1), i.e., there 
are no negative integral exponents in the Cesaro development of /, and consequently 
the constants fj, m are the moments of /. Thus (6.4) simplifies to 



*(t) ~Vc, T(a k + 1) t-"*- 1 + Y" ,^ * 

A — ' *■ — ' ml 

fc>l m>0 



This is precisely the case for a (generalized) Laplacian: if n is odd, only half-integer 
powers of A appear in the spectral density and logarithmic terms in the heat kernel 
are thereby ruled out. Notice that the Cesaro development for an odd dimensional 
Laplacian need not terminate. For even dimensions, the term k = n/2 is proportional 
to wres-ff°A -1 and later terms are proportional to wresi7 r A _r_1 . However, since H r 
is a differential operator, its local Wodzicki residue vanishes for r G N, and the Cesaro 
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development terminates at the A term. However, as we have seen, at this point the 
moments (6.2) enter the picture. 

It has become a habit to write the diagonal of the heat kernel for a Laplacian in 
the form 



K(t, x, x) ~ (47rf)- n / 2 J2 b k(x, x) t k l 



k=0 



where n is the dimension of the manifold and bo(x, x) = 1. We see now that bk(x,x) = 
for k odd, whereas 

h t \ 2 fc a 2fc (x) 

02k(x,x) = -t^—, w r; ; tt tor k > 0. 

v ' ' fi n (n-2)(n-4)...(n-2/c) 

A similar formula holds off-diagonal. As we have noted, these expansions are local in 
the sense that they do not distinguish between a finite and an infinite region of M n , say. 
However, the smallness of the terms after the first is not uniform near the boundary, 
and hence the "partition function" 

K(t):= K(t,x,x)dx~(4nt)- n / 2 J2 b kt k/2 , (6.5) 

Jm k=o 

with bo = vol(M) for scalars, has an expansion with nontrivial boundary terms in 
general, starting to contribute in 62 [33]. 

As for the examples, the expansion (6.5) for § 2 was first obtained as the partition 
function of a diatomic molecule [34] and is well known to physicists. On using vol(§ 2 ) = 
4/T, we read Mulholland's expansion directly by looking at (5.2): 

^W~7 + ^ + ^rt+-^-t 2 + --- as UO. 
t 3 15 315 

As for the SU(2) group manifold, from (5.3), on using vol(S 3 ) = 2-7T 2 and e tA = 
e t e -t(i-A)^ ^e partition function is seen immediately to be 



K&(t) 



4t 3 / 2 



We turn at last to the Chamseddine-Connes expansion. The theory of Cesaro and 
parametric expansions justifies (1.1), in the following way. We work in dimension n = 4 
and take H = D 2 , a generalized Laplacian, acting on a space of sections of a vector 
bundle E, over a manifold without boundary. The kernel of its spectral density satisfies 

rki? 1 

d D 2(x,x;\) ~ ———\ + ——wresD~ 2 (x) (C) as A — > 00. 

lD7T z Szn* 
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Integrating over M and using the formulas of this section with t = A -2 , we then get 

Tr (j>(D 2 /A 2 ) ~ — — - I rk E A 4 / A0(A) rfA + 6 2 A 2 / 0(A) rfA 
(4tt) 2 V. Jo Jo 

+ ]T (-l)"V (m) (0) b 2m+4 (D 2 ) A~ 2 A as U 0. 

m>0 ' 

where (— l) m 6 2 m+4(-D 2 ) = 167r 2 // m (.D 2 )/m! are suitably normalized, integrated moment 
terms of the spectral density of D 2 . Thus, we arrive at (1.1). 

We finally take stock of the status of the Chamseddine-Connes development. If 
(j) G 5, then the development becomes a bona fide asymptotic expansion. However, if 
one wishes to use (for instance) the counting function N^i^X < A 2 ), which does not lie 
in S, then the present formulae are not directly applicable and one one must proceed like 
in Section 5; moreover the expansion beyond the first piece is only valid in the Cesaro 
sense. We close by noting that third piece of the Chamseddine-Connes Lagrangian has 
interesting conformal properties; this is better studied through the corresponding zeta 
function at the origin [35]. That term is definitely not a Wodzicki residue but a moment; 
whether this fact has any physical significance is not easy to say. 
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